Multifractality and scale invariance in human heartbeat dynamics 
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Human heart rate is known to display complex fluctuations. Evidence of multifractality in heart 
rate fluctuations in healthy state has been reported [Ivanov et al., Nature 399, 461 (1999)]. This 
multifractal character could be manifested as a dependence on scale or beat number of the probability 
density functions (PDFs) of the heart rate increments. On the other hand, scale invariance has been 
recently reported in a detrended analysis of healthy heart rate increments [Kiyono et al., Phys. Rev. 
Lett. 93, 178103 (2004)]. In this paper, we resolve this paradox by clarifying that the scale invariance 
reported is actually exhibited by the PDFs of the sum of detrended healthy heartbeat intervals taken 
over different number of beats, and demonstrating that the PDFs of detrended healthy heart rate 
increments are scale dependent. Our work also establishes that this scale invariance is a general 
feature of human heartbeat dynamics, which is shared by heart rate fluctuations in both healthy 
and pathological states. 

PACS numbers: 87.19.Hh,87.10.+e,05.40.-a 



I. INTRODUCTION 

The heartbeat interval in human is known to dis- 
play complex fluctuations, referred to as heart rate 
variability (HRV). In the past decade, many analy- 
ses [3, @, y, EL H, @, 0, Hi have been carried out to char- 
acterize the statistical features of human HRV, with an 
aim to gain understanding of human heartbeat dynam- 
ics. An intriguing finding is the multifractality in healthy 
HRV and the loss of this multifractal character in patho- 
logical HRV in patients with congestive heart failure 
Such multifractal complexity in healthy HRV was fur- 
ther shown to be related to the intrinsic properties of the 
control mechanisms in human heartbeat dynamics and 
is not simply due to changes in external stimulation and 
the degree of physical activity Q. 

In another complicated phenomenon of fluid turbu- 
lence, physical measurements are also known to be mul- 
tifractal [9(. In fluid turbulence, it is common to study 
structure functions, which are the statistical moments of 
the increments of the signals at different scales, and their 
scaling behavior. Multifractality manifests itself as a 
nonlinear dependence of the scaling exponents on the or- 
der of the structure functions. This nonlinear dependence 
is equivalent to the scale dependence of the probability 
density functions (PDFs) of the increments of the signals 
at different scales. These ideas of structure functions in 
fluid turbulence were employed to analyse healthy HRV 
and similar multifractality, with a scale dependence of the 
PDFs of the heart rate increments at different scales or 
between different number of beats, was indeed found 
This analogy of human HRV to fluid turbulence was fur- 
ther exploited and a hierarchical structure found in fluid 
turbulence [h| was shown to exist also in human HRV, 



with different parameters for heart rate fluctuations in 
healthy and pathological states [12] • The different values 
of the parameters can thus be used to quantify the mul- 
tifractal character of healthy HRV and its loss in patho- 
logical HRV. On the other hand, in a recent detrended 
analysis [H| that aims to eliminate the non-stationarity 
of heartbeat data, "scale invariance of the PDFs of de- 
trended healthy heart rate increments" was reported, and 
interpreted as an indication that healthy heartbeat dy- 
namics are in a critical state. This finding appears to be 
in contradiction to the multifractal character of healthy 
HRV and needs clarification. 

In this paper, we resolve this apparent paradox and 
further establish that human heartbeat dynamics exhibit 
a general scale invariance, which is shared by heart rate 
fluctuations in both healthy and pathological states. 

This paper is organized as follows. We first review the 
statistical character of multifractality in healthy HRV in 
Sec. II. In Sec. Ill, we clarify that the scale invariance re- 
ported for healthy HRV in Ref. [l3| is actually exhibited 
by the PDFs of the sum of detrended heartbeat intervals 
and demonstrate explicitly that the PDFs of detrended 
healthy heart rate increments are indeed scale depen- 
dent and thus consistent with the multifractal character 
of healthy HRV. Then we show that pathological heart 
rate fluctuations in patients with congestive heart fail- 
ure also display this scale invariance. Our finding thus 
shows that such scale invariance cannot be an indication 
of healthy human heartbeat dynamics being in a critical 
state. In Sec. IV, we show that this general scale in- 
variance in human heartbeat dynamics is non-trivial in 
that it is absent in multifractal turbulent temperature 
measurements in thermal convective flows. In Sec. V, we 
show that the essential effect of the detrended analysis 
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is to take out the local mean from the data. Finally, we 
summarize and conclude our paper in Sec. VI. 



II. STATISTICAL SIGNATURE OF 
MULTIFRACTALITY 

For completeness, we first review how the multifrac- 
tality of healthy human HRV can be studied using the 
ideas of structure functions in fluid turbulence. Consider 
a dataset of human heartbeat intervals b(i), where i is the 
beat number. The beat-to-beat interval is also known as 
RR interval as it is the time interval between successive 
"R" peaks in the electrocardiogram (ECG) time signal. 
The value of b(i) varies from beat to beat and this varia- 
tion is the human HRV. Following the ideas of structure 
functions in turbulent fluid flows, one defines the 
heart rate increments between an interval of n beats as, 

A n b(i) = b(i + n) - b(i) , (1) 

which arc differences between the heart rate intervals sep- 
arated by n beats. The p-th order structure functions are 
the p-th order statistical moments of the increments: 

S p (n) = <|A n 6(i)|*) (2) 

It was found [l(| that S p (n) for heart rate fluctuations 
in healthy state exhibits power-law dependence on n: 

S p (n) ~ nfc (3) 

for n w 8 — 2048. This power-law or scaling behavior 
is analogous to that found for velocity or temperature 
structure functions in turbulent fluid flows. Moreover, 
the scaling exponents Q p were found to depend on p in a 
nonlinear fashion [101 ] as in turbulent fluid flows. This im- 
plies that for healthy HRV, the standardized PDFs (with 
mean substracted then normalized by the standard de- 
viation) of Ab n changes with the scale or the number 
of beats n, and are thus scale dependent. In fact when 
Eq. © holds, the standardized PDFs of Ab n are scale 
invariant, i.e., independent of n, if and only if Q p is pro- 
portional to p. Hence the nonlinear dependence of £ p on 
p or, equivalently, the scale dependence of the standard- 
ized PDFs of Ab n on n is a characteristic signature of the 
multifractality of healthy HRV, in analogy to turbulent 
fluid flows. 



III. SCALE INVARIANCE IN THE 
DETRENDED ANALYSIS 

In contrast to physical measurements in turbulent fluid 
flows, human heartbeat interval data are often non- 
stationary. This non-stationarity is one possible reason 
for the relatively poor quality of scaling in HRV as com- 
pared to that in turbulent fluid flows. To eliminate the 
non-stationarity, a "detrended fluctuation analysis" has 



been introduced [14]]. This analysis was further devel- 
oped to study detrended heart rate increments [H, [H[ . 
The procedure of this detrended analysis consists of the 
following steps. First, B(m), which is the sum of b(j): 

m 

B(m) = X>(j) , (4) 
i=i 

is calculated. Second, the dataset of B(m) is divided into 
segments of size 2n, and the datapoints in each segment is 
fitted by the best gth-order polynomial. This polynomial 
fit represents the "trend" in the corresponding segment. 
Third, these polynomial fits are subtracted from B{m) 
to get B* (m) , which are then "detrended" . Finally, the 
standardized PDFs of the increments of B*: 

A n B*(i) = B*(i + n) - B*{i) (5) 

for different values of n are studied. Note that in 
Refs. [H, HU, A n B*(i) was denoted as A n B(i) and the 
symbol s was used in place of n. The standardized PDFs 
of A n B*(i) for healthy heartbeat data were found to 
be independent of n, and this was referred to as "scale- 
invariance in the PDFs of detrended healthy human heart 
rate increments" in Ref. [l3| ■ We shall show below that 
this conclusion is inaccurate. 

Let us denote the detrended heartbeat interval by 
b*(i). From the detrended procedure described above, 
it is natural to define b*(i) by 

m 

B*(m)=J2b*(j) (6) 
i=i 

Then 

i+n 

A n B*(i) = ]T b*(j) (7) 
j=i+i 

Thus the detrended heart rate increment between n beats 
should be defined as 

A n b*(i) = b*(i + n)-b*(i) (8) 

Hence A n B*(i) is the sum of detrended heartbeat in- 
tervals taken over n beats rather than detrended heart 
rate increments. As a result, the observation of scale- 
invariant or n- independent standardized PDFs of A n B* 
does not necessarily imply that the standardized PDFs 
of A n b* are also n-independent. Indeed, one expects the 
contrary, namely, the standardized PDFs of A n b* should 
depend on n as healthy human HRV is multifractal. 

To investigate this issue, we study the scaling behav- 
ior of the statistical moments of A n B* and A n b* . As 
the standardized PDFs of A n B* are n independent, the 
scaling exponents for the statistical moments of A n B* 
should be proportonal to the order of the moments. On 
the other hand, we expect the scaling exponents for the 
statistical moments of A n b* to have a nonlinear depen- 
dence on the order of the moments. We analyze healthy 
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heartbeat data that are taken from a database of 18 sets 
of daytime normal sinus rhythm data downloaded from 
public domain[I3|. We follow the detrended procedure 
described above to get A n B*(i). We find that a polyno- 
mial of degree 3 (q = 3) is sufficient to fit the "trend" 
as found in Ref. [la ]. To get the detrended heart rate 
increment b*(i), we use Eq. ^ to get 



b*(i) = B*(i) — B*(i — 1) 



(9) 



for both B*(i — 1) and B*(i) belonging to the same seg- 
ment and skip that datapoint when B*(i — 1) and B*(i) 
fall into different (consecutive) segments. Next, we eval- 
uate the statistical moments 



S p (n) 
S*(n) 



(|A„B*(0n 

(|a„&*«h 



(10) 

(11) 



As seen from Fig. [TJ S p (n) exhibits power-law or scaling 
behavior with n with exponents C p : 



Sp(n) 



(12) 



for n between 16 to 1024 and p between 0.2 to 3. On the 
other hand, S*(n) exhibits better scaling behavior with 
n with exponents C*: 



s;(n) 



(13) 



for n between 32 to 1024 and p between 0.2 to 3 (see 
Fig. [I]). 

In Fig. [31 we plot the scaling exponents ( p and Cp as a 

function of p. It can be seen that Cp is proportional to p 
confirming that the standarized PDFs of A„ B* are scale 
invariant as reported in Ref.[l3J]. On the other hand, Cp 
is not proportional to p but changes with p in a nonlin- 
ear manner, as expected from the multifractal character 
of healthy human HRV. Thus we have clarified that for 
healthy human HRV that is multifractal, the standard- 
ized PDFs of detrended heart rate increments between n 
beats depend on n while those of the sum of detrended 
heartbeat intervals taken over n beats are n-independcnt. 
We have also calculated the scaling exponents C p of S p (n) 
of the statistical moments of untreated heart rate incre- 
ments and the results are shown in Fig. [3] too. It can be 
seen that the detrended procedure does not change much 
the scaling exponents of the heart rate increments. We 
shall return to this in Sec. V. 

It was suggested that this scale invariance of the stan- 
dardized PDFs of A n B*, the sum of detrended heartbeat 
intervals taken over n beats, is an indication of hea lthy 
human heartbeat dynamics being in a critical state [1 31 ] - 
To check this suggestion, it would be useful to perform 
the same analysis to human HRV in pathological state. 
Thus, we perform the same analysis using 45 sets of 
daytime data from congestive heart failure patients, also 
downloaded from the same public domain [l6j . The re- 
sults for C p and C* in this case are shown in Fig. [|J Note 
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FIG. 1: The statistical moments S p (n) of the sum of de- 
trended heartbeat intervals [see Eq. (I10|l for definition] for 
healthy heartbeat data for p — 0.2 (circles), p = 0.6 (trian- 
gles), p = 1.0 (squares), p = 1.6 (plusses), p = 2.0 (crosses), 
p — 2.6 (diamonds), and p — 3.0 (inverted triangles). The 
curves have been shifted vertically for clarity. 



that C* is now approximately proportional to p, show- 
ing that the multifractality is lost in pathological HRV, 
consistent with earlier report On the other hand, C P 
is again proportional to p, demonstrating that the scale 
invariance of the standardized PDFs of the sum of de- 
trended heartbeat intervals is not restricted to healthy 
HRV but also exhibited by pathological HRV in conges- 
tive heart failure patients. Moreover, we find that the 
scale- invariant standardized PDFs are approximately ex- 
ponential for both the healthy and pathological heartbeat 
data as shown in Figs. O and O 

Since this scale invariance is found generally in heart 
rate fluctuations in both healthy and pathologial state, 
it could not be an indication that healthy heartbeat dy- 
namics are in a critical state. Common feature for both 
healthy and diseased human HRV was also reported be- 
fore [20| : it would be interesting to explore whether this 
earlier feature and the present one are related. 



IV. DETRENDED ANALYSIS FOR 
TURBULENT TEMPERATURE 
MEASUREMENTS 

It is natural to ask whether this general scale invari- 
ance found in human heartbeat dynamics is trivial, i.e., 
whether it exists for any fluctuating data. In this sec- 
tion, we shall see that such scale invariance is absent in 
temperature data in turbulent flows so the answer to the 
above question is no. 
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FIG. 2: The statistical moments Sp(n) of detrended heart 
rate intervals [see Eq. |TT} for definition] for healthy heartbeat 
data for p ranges from 0.2 to 3.0. Same symbols as in Fig. 1. 
The curves have been shifted vertically for clarity. 



FIG. 4: The exponents Cp/C2 (plusses) and Cp/CJ (crosses) as 
a function of p for pathological heartbeat data from congestive 
heart failure patient. Both of them are close to the solid line 
of p/2. 




FIG. 3: The exponents Cp/C2 (plusses), Cp/CI (crosses), and 
Cp/Ca (circles) as a function of p for healthy heartbeat data. 
It can be seen that Cp/& is close to p/2 which is shown as the 
solid line. 



Specifically, we apply the detrended analysis to tem- 
perature measurements taken in turbulent thermal con- 
vective flows [ljj]- In place of b(i), we now have 0(ti), the 
temperature measurement taken at time ij. In the ex- 
periment, the measurements were sampled at a constant 
frequency of 320 Hz such that ti = iSt with 8t = 1/320 s. 



The standardized PDFs of the temperature increments 
A n 0(ti) = 9{tj+n) — 0(ti) have been studied and found to 
change with n [1 81 ] thus the temperature data in turbulent 
thermal convection are multifractal. Also, the tempera- 
ture structure functions, R p (n) = (\A n 8(ti)\ p ) have been 
studied and found to have good relative scaling (l9j : 

R p (n) ~ [i? 2 (n)]^/« 2 (14) 

We calculate 0(i m ) = Y^jLi ^fe) an d repeat the pro- 
cedure of the detrended analysis, as discussed in Sec. Ill 
with B(m) replaced by 0(i m ), to obtain A n Q*(ti) and 
A n 0*{ti). We then calculate the corresponding statis- 
tical moments R p (n) = (|A„6*(^)| P ) and R*(n) = 
(\A n 6* (ti)\ p ) and their respective relative exponents 
£ p /£ 2 and defined by: 

Rp(n) ~ [R 2 (n)f^ (15) 
R* p (n) ~ [JJ3(n)]Q/« (16) 

Our results are shown in Fig. [7J Again we find that 
£p/£! deviates from p/2, as expected from the multi- 
fractal character of the turbulent temperature measure- 
ments. However, interestingly £ p /£;2 deviates from p/2 
too, showing that the standardized PDFs of A„8* are 
scale dependent and changing with n. To show this de- 
viation more clearly, we plot — p/2 versus p in the 
inset of Fig. El 

Indeed, the standardized PDFs of A n O* changes from 
stretched-exponential to exponential to Gaussian as n in- 
creases from 4 to 4096, as shown explicitly in Fig.[5J This 
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FIG. 5: Standardized PDFs for A n B* for healthy heartbeat 
data with n = 4 (circles), n = 16 (squares), n = 64 (dia- 
monds), and n = 256 (triangles). Data from four different 
healthy subjects are shown and seen to coincide with one an- 
other. These n-independent PDFs are seen to be well approx- 
imated by the standard exponential distribution (solid line). 



change of the standardized PDFs of A„9*, the sum of 
detrended temperature measurements taken over n sam- 
pling intervals, with n is similar to the change of the 
standardized PDFs of the temperature increments A n 6 
with n as reported in Ref. [HI]. 

As can be seen in Fig. are c l° se to £ P /£2, indi- 

cating again that the detrended procedure does not affect 
the scaling exponents of the temperature increments. We 
shall understand this in the next section. 



V. THE ESSENTIAL EFFECT OF THE 
DETRENDED ANALYSIS 

In this section, we shall explore and understand what 
the detrended procedure does to the data. As discussed 
in Sec. Ill, the "trend" in the data is estimated by a 
polynomial fit in each segment of the dataset of B(m), 
and we have used a polynomial of degree 3. We check 
that our results do not change much when a polynomial 
of a lower degree is used instead. In particular, we ob- 
tain similar results by using a linear fit of the different 
segments of B(m). In the following, we shall get explicit 
results for "detrended" B* when the "trend" is estimated 
by a linear fit. 

Let us focus on the Zth segment of B(m) with m\ < 
m < TO2, where mi = (I — l)(2n) + 1 and rri2 — l(2n) 
for some I. I runs from 1, 2, 3, . . . for all the segments. 
Denote the best linear fit to this segment by a; m + c/ 
where the fitting constants a; and c/ depend on I. The 



FIG. 6: Standardized PDFs for A n B* for heartbeat data from 
a congestive heart failure patient with n — 4 (circles), n 
I f! (squares), n — 64 (diamonds), and n = 256 (triangles). 
Again, the scale invariant PDFs are well approximated by 
the standard exponential distribution (solid line). 



fitting constant a; can be reasonably well approximated 
by the slope in this segment: 



B(m 2 ) - g(mi) 
2n- 1 



Using Eq. (g]), we have 



ai 



2n - 1 



2n 



(17) 



k (18) 



where bi is the local average of b(j) in the Ith segment. 
Recall from Sec. Ill that B* is B subtracting the best 
linear fit and use Eq. ([4]), we have 



B*(m)«5>C7)-5l]-Ci 



(19) 



and 



B* (m + n) - \ Hj) ~ 1{ \- Cl m + n<m 2 

a \m -h n) ~ < m+ „ 

I L,j=i [b{J) - bi+i\ - q+i m + n > m 2 



Thus using Eq. (O, we get 



j=m+l 



(20) 
(21) 



for m + n < m,2 and 

m,2 



A„B*(m)« Y, \PU)-bl]+ E P»C3)-&l+i] ( 22 ) 

j=m+l j=m 2 + l 
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FIG. 7: The three relative exponents (, P /(,2 (plusses), 
ip/Q (crosses), and £ P /£2 (circles) for temperature measure- 
ments in turbulent convective flows. All the three relative 
exponents deviate from p/2 (the solid line). The deviation 
t,p/t,2 — p/2 is plotted versus p in the inset to show clearly 
that £p/£2 is not proportional to p. 

for m + n > m^- To obtain Eq. (f2"2")l . we make use of 
the approximation that the two linear fits of the /th and 
(/ + l)th segments intersect at m = ni2- 

bi m 2 + ci w b t+ i 772 2 + cj + i 
=> Q + i - Q » (6j - bi +1 )m 2 (23) 

Comparing Eqs. (|21[) and (f2"2"|) with ([7]), we see immedi- 
ately that the detrended heartbeat interval b* is given 
approximately by 

( 24 ) 

Hence what the detrended analysis essentially does is to 
subtract the local average from the data. 

To verify this directly, we redo the analysis for the 
heartbeat data with the local mean subtracted and com- 
pare the results obtained with those from the detrended 
analysis. We define 

b(j) = b{j)-bi (25) 

m 

B(m) = J2Hi) (26) 

8=1 

and study the scaling behavior of the statistical moments 
of 

A„S(j) = b(j + n)-b(j) (27) 

j+n 

A n B(j) = B(j + 77) - B(j) = J2 ~ h V> (28) 

i=j+i 




-6 -4 -2 2 4 6 
A 0* 

n 

FIG. 8: The standardized PDFs of A n 0* for temperature 
measurements taken in turbulent thermal convective flows 
with n = 4 (solid), n = 32 (dashed), n = 256 (dot-dashed) 
and n = 4096 (dotted). The dependence of the standardized 
PDFs on n is clearly seen. 

with n. The corresponding exponents are denoted by a p 
and dp, which are defined by 

(\A n b(i)\P) ~ n a - (29) 
<|A n B(i)|*) ~ n*> (30) 

We compare a p and a p with and Cp respectively. As 
seen from Fig. [5J the results are in good agreement con- 
firming that the essential effect of the detrended analysis 
is to eliminate the local average from the data. 

As a result, b*(j + n) — b*(j) w b(j + n) - b(j) will 
be close to b(j + n) — b(j), that is, the heart rate incre- 
ments are not affected much by the detrended analysis. 
This explains why (* are close to ( p (see Fig. [3]) and 
similarly why £* are close to £ p (see Fig. [7]). On the 

other hand, B* (j + n) - B* (j) w B(j + n) - B(j) can be 
different from B(j + n) — B(j), and thus the sum of de- 
trended heartbeat intervals could have different statisti- 
cal features from those of the sum of untreated heartbeat 
intervals. 



VI. SUMMARY AND CONCLUSIONS 

Understanding the nature of the complicated human 
HRV and thus human heartbeat dynamics has been the 
subject of many studies. An interesting and intrigu- 
ing finding reported in earlier studies [s| is that in the 
healthy state, human heart rate fluctuations display mul- 
tifractality, and that this multifractal character is lost 
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FIG. 9: Comparison of the exponents Cp/C2 (plusses) and 
CpKZ (crosses) with a p /&2 (circles) and a v /a.2 (squares) ob- 
tained respectively from the detrended analysis and from the 
analysis eliminating the local mean. The comparison for 
healthy heartbeat data is shown in (a) and (b) while that 
for pathological heartbeat data from congestive heart failure 
patients is shown in (c) and (d). Good agreement between ( p 
and & p and between ^ and a p is seen. 



intervals at different scales or between different number 
of beats. A detrending analysis aiming to eliminate the 
non-stationarity of heartbeat data has been performed, 
and " scale invariance of the PDFs of detrended healthy 
human heart rate increments" reported [l3j]. We have 
resolved this paradox by clarifying that the scale invari- 
ance found in the detrended analysis is actually exhibited 
by the PDFs of the sum of detrended heartbeat intervals 
taken over different number of beats, and demonstrat- 
ing explicitly that the PDFs of detrended healthy heart 
rate increments are scale dependent. We have under- 
stood the essential effect of this detrended analysis is to 
eliminate the local average from the heartbeat data. We 
have further found that this scale invariance of the PDFs 
of the sum of heartbeat intervals, with the local mean 
subtracted, is displayed also by heart rate fluctuations 
of congestive heart failure patients. In both the healthy 
and pathological states, such scale-invariant PDFs are 
close to an exponential distribution. On the other hand, 
this scale invariance is absent in the multifractal tem- 
perature measurements in turbulent thermal convective 
flows. Hence we have found an interesting scale invari- 
ance of exponential PDFs in human heartbeat dynamics, 
which is exhibited generally by heart rate fluctuations in 
both healthy and pathological states. Since this scale in- 
variance is a general feature, it cannot be an indication 
of the healthy state being critical, in contrast to what 
was claimed in earlier studies 1131. 



for heart rate fluctuations in pathological state such as 
congested heart failure. Based on an analogy with mea- 
surements in turbulent fluid flows, which are known to 
have multifractal character as well, such multifractality 
in healthy HRV can be manifested as a scale dependence 
of the standardized PDFs of the increment of heartbeat 
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